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Abstract — The  self-consistent  time-dependent  code  MAGY  is 
presently  being  used  for  designing  and  modeling  gyrodevices. 
In  this  paper,  the  code’s  self-consistent  capability  is  used  to 
investigate  three  different  issues  relating  to  the  operation  of 
gyroklystron  amplifiers.  These  are  the  effect  of  window  reflection 
on  the  properties  of  the  output  waves,  higher  order  mode  exci¬ 
tation  in  nonlinear  output  tapers,  and  excitation  in  cutoff  drift 
sections.  The  first  two  effects  can  potentially  have  major  impacts 
on  the  gyroklystron  output  radio  frequency  (RF)  in  terms  of 
bandwidth  power  ripple,  phase  ripple  across  the  band,  and  mode 
purity.  In  fact,  strong  correlation  between  simulation  ^results  and 
experimental  data  indicates  that  the  observed  bandwidth  ripple 
in  a  recently  tested  W-band  gyroklystron  can  be  attributed  to  a 
small  (— 20-dB)  window  reflection.  The  last  effect,  excitation  in 
cutoff  drift  sections,  is  found  to  result  in  a  high  level  of  RF  fields  at 
the  drive  frequency  in  drift  sections.  This  effect  impacts  the  design 
and  use  of  lossy  ceramics  in  the  drift  section  because  of  thermal 
considerations,  particularly  in  high  average  power  devices. 

Index  Terms — Gyroamplifiers,  gyroklystrons,  output  tapers, 
window  reflection. 


I.  Introduction 

DETAILED  and  accurate  nonlinear  simulations  of  the 
interaction  between  the  electromagnetic  (EM)  fields  and 
the  electron  beam  in  vacuum  electron  devices  are  essential 
to  predict  their  performance.  The  prediction  capabilities  are 
critical  to  the  design  and  development  of  future  devices. 
This  is  the  essential  impetus  for  the  recent  development  of 
the  self-consistent,  time-dependent,  quasi  three-dimensional 
(3-D),  new  MAGY  code  [1].  In  contrast  to  the  finite-differ¬ 
ence-time-domain  particle-in-cell  (PIC)  approach,  MAGY 
employs  a  reduced  description  approach  in  which  the  EM  fields 
are  described  by  a  superposition  of  the  waveguide  transverse 
electric  (TE)  and  transverse  magnetic  (TM)  eigenmodes. 
Furthermore,  the  temporal  evolution  of  the  EM  complex  field 
amplitudes  and  of  the  electron  beam  are  assumed  to  be  slow 
relative  to  the  radio  frequency  (RF)  period,  hence,  allowing  for 
the  averaging  over  time  scales  on  the  order  of  the  RF  period. 
The  combination  of  fast  time-scale  averaging  and  the  reduced 
description  RF  fields  substantially  reduce  the  computational 
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resources  compared  with  that  required  by  PIC  codes.  However, 
unlike  other  reduced  description  codes,  MAGY  does  not 
incorporate  restricting  assumptions  on  the  physics  involved 
and,  thus,  does  not  compromise  on  the  fidelity  of  the  results. 
For  instance,  in  MAGY,  the  RF  field  profiles,  rather  than  being 
restricted  to  a  fixed  form,  freely  evolve  in  response  to  the 
interaction  with  the  electron  beam.  This  capability  is  achieved 
via  a  novel  formulation  of  the  generalized  telegrapher’s  equa¬ 
tions  together  with  the  fast  time-scale  averaging  of  the  fully 
relativistic  electron  equations  of  motion,  which  provides  for 
multimode  (TE  and  TM)  coupling  in  arbitrary  wall  radius  pro¬ 
files,  especially  at  radial  step  discontinuities.  This  capability, 
in  single  mode  form,  has  previously  been'discussed  in  [2].  The 
new  MAGY  code  is  intended  for  the  modeling  of  both  slow- 
and  fast- wave  devices.  However,  in  this  paper,  we  shall  only 
report  results  on  a  type  of  fast- wave  devices — gyroklystron 
amplifiers. 

Specifically,  MAGY  is  employed  to  investigate  three  dif¬ 
ferent  effects  in  gyroklystrons.  First  is  the  effect  of  window 
reflection  on  the  properties  of  the  output  waves.  The  excellent 
agreement  between  experimental  data  and  MAGY  simulations 
strongly  suggests  that  the  observed  bandwidth  power  ripple 
from  the  recently  completed  test  of  a  W-band  gyroklystron  [3], 
[4]  is  caused  by  small  power  reflections  (—20  dB)  from  the 
output  window.  In  addition  to  the  power  ripple,  it  can  also  be 
shown  that  reflections  introduce  ripple  on  the  output  phase  and 
modify  the  effective  quality  factor  of  the  output  cavity  across 
the  band.  This  effect  is  discussed  in  Section  n.  The  excitations 
of  higher  order  modes  in  the  output  taper  are  presented  in 
Section  III.  The  internal  mode  conversion  process  inside  the 
output  taper  provides  the  initial  seed  for  power  in  higher  order 
modes,  which  are  then  amplified  via  beam-wave  interactions  of 
the  gyro  traveling- wave-tube  type  (gyro-TWT)  in  the  uptaper. 
These  interactions  can  compromise  the  high  mode  purity  of  the 
output  waves  for  which  the  taper  is  designed.  Third,  Section  IV 
discusses  the  existence  of  RF  excitation  at  the  drive  frequency 
in  cutoff  drift  sections  driven  by  the  bunched  electron  beam. 
This  phenomenon  may  impact  the  design  and  use  of  lossy 
ceramics  in  drift  sections  of  high-average  power  devices 
because  of  thermal  considerations.  Section  V  summarizes  and 
concludes. 

II.  Effect  of  Output  Window  Reflection 

It  is  well  established  that  reflection  of  microwaves  from  the 
output  window  can  influence  the  performance  of  gyrotron  os¬ 
cillators  [5]-[9].  Small  window  reflections  can  result  in  mode 
competition  and  frequency  jumping  in  these  devices.  In  gyroam¬ 
plifiers,  window  reflections  are  responsible  for  the  power  ripple 
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Fig.  1.  Experimental  bandwidth  plots  of  the  high  average  power  W-band  gyro-klystron  at  four  different  magnetic  field  settings.  Beam  current  and  voltage  are  4 
A  and  65  kV,  respectively. 

across  the  bandwidth  of  the  high-average  power  W-band  gy- 
roklystron  amplifier  [3],  [4].  The  experimentally  observed  band¬ 
width  ripple  is  illustrated  in  Fig.  1  for  four  different  circuit  mag¬ 
netic  fields.  The  beam  current  and  voltage  for  the  cases  shown 
are  4  A  and  65  kV,  respectively.  As  shown  in  Fig.  1 ,  the  peaks 
and  valleys  of  the  ripple  occur  at  constant  frequencies  for  all 
four  magnetic  fields.  Similar  behavior  has  also  been  observed 
at  all  beam  currents.  The  modulation  of  the  gyroklystron  fre¬ 
quency  response  seen  in  Fig.  1  is  indeed  surprising,  especially 
considering  that  the  reflected  power  is  at  least  20  dB  down  from 
the  outgoing  RF  waves  according  to  window  cold  tests.  In  this 
section,  we  shall  first  discuss  a  simple  analytical  model,  which 
describes  the  effect  of  window  reflection  on  gyroklystron  ampli¬ 
fier  performance.  The  resulting  scaling  laws  are  then  compared 
with  MAGY  simulation  results.  This  section  concludes  with  a 
detailed  comparison  of  MAGY  results  with  experimental  data. 

The  Naval  Research  Laboratory  (NRL)/Litton/Communica- 
tion  and  Power  Industries’  (CPI)  high-average  power  W-band 
gyroklystron  amplifier  is  a  stagger-tuned  four-cavity  TE01  gy¬ 
roklystron  operating  at  the  fundamental  cyclotron  harmonic  [3], 

[4],  The  amplifier  geometry  is  illustrated  by  Fig.  2,  which  shows 
the  four  cavities  and  the  nonlinear  output  uptaper  (~10.1  cm  in 
length).  The  cavities  are  isolated  from  each  other  by  drift  sec¬ 
tions,  which  are  cutoff  to  the  operating  mode.  Following  the  up¬ 
taper  is  the  collector,  which  also  serves  as  the  RF  output  wave¬ 
guide  [10].  At  the  end  of  the  collector  (~63.6  cm  downstream 
from  the  end  of  the  uptaper)  is  a  single-disk  chemical  vapor 
deposition  (CVD)  diamond  window.  As  alluded  to  earlier,  this 
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Fig.  2.  Geometry  and  a  typical  RF  field  profile  from  MAGY  simulations  of 
the  high  average  power  W-band  gyroklystron.  A  window  reflection  of  —  20-dB 
has  been  assumed. 

window  has  been  cold-tested  and  found  to  have  a  return  loss 
of  approximately  —20-dB  reflection  (1%  power  reflection)  over 
the  frequency  range  of  interest.  Also  shown  in  Fig.  2  is  a  typical 
RF  field  profile  simulated  with  MAGY  with  a  -20-dB  window 
reflection.  This  field  profile  has  two  particular  noteworthy  fea¬ 
tures.  One  is  the  standing  wave  pattern  seen  in  the  nonlinear 
uptaper,  which  results  from  the  interference  between  the  for¬ 
ward  going  wave  and  the  backward-reflected  RF.  The  other  is 
the  presence  of  finite  RF  amplitude  in  the  cutoff  drift  section 
between  the  penultimate  and  the  output  cavities.  This  shall  be 
the  topic  of  discussion  in  Section  IV. 

In  gyroklystrons,  an  annular  gyrating  electron  beam,  usually 
produced  with  a  magnetron  injection  gun  (MIG),  is  introduced 
into  the  interaction  circuit,  which  comprises  of  at  least  two  cav¬ 
ities.  At  the  input  cavity,  the  electron  energies  are  modulated  by 
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Fig.  3.  (a)  Output  power  versus  frequency  with  and  without  window  reflection  and  (b)  corresponding  power  ripple,  A P/P0,  for  constant  RF  current  source 
excitation  of  output  cavity. 


the  TE  fields  introduced  by  the  RF  input  power,  changing  their 
relativistic  cyclotron  frequencies.  The  modulation  of  relativistic 
cyclotron  frequencies  causes  the  electrons  to  bunch  in  phase 
angle.  The  azimuthal  bunching  process  is  further  enhanced  via 
ballistic  bunching  in  the  drift  sections  and  by  the  intermediate 
cavities.  Finally,  the  RF  power  is  diffractively  coupled  from  the 
output  cavity  and  flows  toward  the  output  window.  In  steady 
state,  energy  balance  dictates  that  power  extracted  from  the  elec¬ 
trons  must  equal  the  RF  power  diffractively  coupled  from  the 
cavity  plus  the  ohmic  power  lost  inside  the  cavity.  Assuming 
negligible  ohmic  loss,  the  above  statement  can  be  expressed  as 

PoM  =  ^-W(u)  =  [  ER F(w,f)  •  Jb(w,r)  dV.  (1) 
Vdiff  Jv 


Here,  <x>  =  2irf  is  the  angular  frequency  of  the  input  power. 
Qdiff  is  the  diffractive  quality  factor  of  the  output  cavity.  Erf 
is  the  cavity  RF  electric  field.  Jb  is  the  electron  bunching  current 
density  at  frequency  t o  and  a  function  of  the  RF  electric  field, 
and  W  is  the  output  cavity  stored  energy.  The  stored  energy  is 
given  by 

W(w)  =  |  Erf( w,  ?)  ■  E*rf(w,7)  dV  (2) 

where  the  integral  is  over  the  output  cavity  volume. 

It  is  worthwhile  to  note  at  this  juncture  that  the  output  window 
reflection  will  only  influence  the  output  cavity  because  the  cavi¬ 
ties  are  isolated  from  one  another  by  the  drift  sections.  Defining 
Erf  =  Eq  as  the  output  cavity  RF  electric  field  without  window 
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reflection,  the  cavity  electric  field  with  finite  reflection  coeffi¬ 
cient,  p,  can  be  expressed  as 

Erf  =  E0(1  +  R){1+R  +  R2  +  ---  +  Rn}  (3) 

where  R  =  pc'°  is  the  complex  reflection  coefficient  and  6  is 
the  total  phase  shift  for  the  round-trip  time  from  the  cavity  to 
the  window  and  back.  This  phase  shift  is  given  by  [6] 

L)  =  4>r  ±  2  f  kz(u),z)  dz  (4) 

Jo 

where  <f>r  is  the  phase  shift  induced  by  the  window,  kz  is  the  axial 
wavenumber,  and  L  is  the  distance  between  the  output  cavity 
and  the  window  (including  the  nonlinear  output  taper). 

In  (3),  the  geometric  series  represents  the  contribution  from 
previous  reflections.  Because  p  <  1,  this  series  quickly  con¬ 
verges  and  we  obtain 

=  (5) 

In  the  event  that  Ji  is  weakly  dependent  or  independent  of  Erf 
(e.g.,  linear  regime),  it  follows  from  (1)  that 


Fout(W,L)SP0(w,L)Re 

=  jP°  (l  +  p2  Zjpeos^)  (6) 

where  Po  is  the  output  power  without  reflection.  Because  ^  is  a 
function  of  cu>,  a  bandwidth  power  ripple  on  the  order  of 


can  be  expected.  Using  the  waveguide  dispersion  relation  for 
the  region  between  the  output  cavity  and  window,  the  ripple 
frequency  spacing  is  found  from  (4)  to  be 


A/(w,  L)  = 


J 


2lo  dz 
kz(u,z) 


r  2  dz 

Jo 


(8) 


*) 


where  vg  is  the  group  velocity.  Thus,  for  example,  we  can  expect 
a  bandwidth  power  ripple  of  approximately  ±20%  for  a  — 20-dB 
(p  =  0.1)  reflection.  As  previously  mentioned,  (7)  is  valid  only 
when  Ji  is  weakly  dependent  on  or  independent  of  £rf,  such 
as  near  the  band  edges  or  in  the  linear  regime  where  the  ef¬ 
ficiency  is  low.  At(or  near  saturation,  strong  interaction  with 
the  RF  fields  rapidly  depletes  the  electron-beam  transverse  ki¬ 
netic  energy.  Consequently,  variation  of  the  RF  fields  will  have 
a  weaker  impact  on  the  output  power  level  (well  known  in  am¬ 
plifiers  as  gain  suppression);  hence,  (7)  only  serves  as  an  upper 
limit  for  the  power  ripple. 

To  study  the  impact  of  window  reflection  on  gyroamplifiers 
self-consistently,  a  window  reflection  model  has  been  imple¬ 
mented  in  MAGY.  In  order  to  test  this  model  and  to  validate 
(7)  and  (8),  MAGY  simulations  are  first  performed  by  exciting 
the  output  cavity  with  a  constant  RF  current  source  (without  an 
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Fig.  4.  MAGY  output  power  versus  reflection  coefficient  (dots)  and  expected 
output  power  according  to  (6)  at  the  fixed  frequency  bf  93.595  GHz  for  constant 
current  excitation  of  output  cavity. 

electron  beam).  The  source  has  a  Gaussian  axial  profile.  This  is 
the  technique  to  determine  the  cold-cavity  resonance  frequency 
and  cavity  quality  factor  (Q)  in  MAGY.  Shown  in  Fig.  3(a)  are 
the  output  powers  versus  frequency  with  (-20  dB)  and  without 
window  reflection.  The  power  curve  without  reflection  shows 
the  output  cavity  cold  resonance  frequency  to  be  93.62  GHz 
with  a  —  3-dB  bandwidth  of  0.567  GHz  corresponding  to  Qaiff 
=  165.  More  illustrative  is  Fig.  3(b),  which  shows  the  power 
ripple  (AP/Pq)  versus  frequency  for  the  two  power  curves  in 
Fig.  3(a).  The  ripple  amplitude  is  approximately  ±20%,  as  ex¬ 
pected  from  (7).  The  ripple  spacing,  A/,  is  0.17  GHz  and  is 
consistent  with  (8)  for  L  &  73.6  cm  (including  the  nonlinear 
uptaper  shown  in  Fig.  1).  The  scaling  of  output  power  versus  re¬ 
flection  coefficient,  /?,  at  a  fixed  frequency  is  illustrated  in  Fig.  4 
for  both  MAGY  simulations  and  (7)  with  excellent  agreement. 
The  frequency  for  this  particular  case  is  93.595  GHz,  which  cor¬ 
responds  to  <j)  =  0,  as  can  be  inferred  from  Fig.  3.  We  note  here 
that  MAGY  simulations  do  not  model  the  device  all  the  way 
to  the  output  window,  which  is  computationally  intensive.  It  is 
only  necessary  to  induce  a  reflection  at  the  end  of  the  simulation, 
which  properly  accounts  for  the  window  reflection  coefficient, 
the  delay  time,  and  the  phase  shift  introduced  by  the  waveguide 
in  accordance  with  (4).  For  the  case  at  hand,  the  end  of  the  sim¬ 
ulation  is  also  the  end  of  the  uptaper,  as  illustrated  in  Fig.  2. 

Recall  from  (1 )  and  (2)  that  although  the  output  power  is  ap¬ 
proximately  proportional  to  2?rf,  the  cavity  stored  energy  is 
proportional  to  |£rf|2-  Thus,  the  cavity  stored  energy  with  re¬ 
flection  is  given  by  the  stored  energy  without  reflections,  Wo, 
by 

<9) 

This  scaling  is  borne  out  by  Fig.  5,  which  shows  the  cavity  stored 
energy  with  and  without  reflection  corresponding  to  the  cases 
shown  in  Fig.  3.  The  main  implication  of  the  scaling  given  by 
(7)  and  (9)  when  taken  together  is  the  modification  of  the  output 
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Fig.  5.  Output  cavity  stored  energy  versus  frequency  corresponding  to  the  constant  current  excitation  case  shown  in  Fig.  3. 


cavity  diffractive  quality  factor  .  Consequently,  in  the  pres¬ 
ence  of  reflection  the  effective  cavity  quality  factor,  Qeff,  is 
given  by 


Qeff(w,L)  = 


uW{u,L) 

1  +  p2  +  2p  cos  <; 


=  Qdifr 


i -P2 


(10) 


which  is  a  function  of  both  frequency  and  reflector  position. 
This  implies  that  small  reflections  can  significantly  impact  the 
start  oscillation  current  threshold  of  the  output  cavity  in  gy- 
roklystrons. 

In  addition  to  causing  bandwidth  power  ripple,  window  re¬ 
flections  also  induce  a  ripple  in  the  output  phase  across  the  band 
relative  to  the  output  phase  without  reflections.  It  can  be  ex¬ 
pected  from  the  above  discussions  that  the  relative  output  phase 
ripple  induced  by  reflection,  A <p,  will  be  identically  zero  at 
the  peaks  and  valleys  of  the  bandwidth  power  ripple  (where 
</>  =  n7r,  n  =  0, 1, 2,  3  •  *  *  etc.).  The  relative  output  phase  ripple 
will  be  at  a  maximum  or  a  minimum  when  the  round-trip  phase 
shift  is  (f>  =  mr/2,  n  =  1,  3, 5  •  •  *  etc.  At  these  points,  assuming 
the  current  source  is  independent  of  Erf,  such  as  in  the  constant 
current  source  excitation,  the  reflected  wave  does  not  interact 
with  the  current  source  in  the  cavity  and  is  simply  reflected  back 
to  the  output  waveguide.  From  power  balance  considerations, 
the  relative  output  phase  is 


A<p(oj,  L )  =  Atan  ((1^^1/2)  •  0  D 


The  output  phase  versus  frequency  corresponding  to  the  cases 
shown  in  Fig.  3  from  MAGY  simulations  is  plotted  in  Fig.  6(a). 
The  difference  between  the  output  phase  with  and  without 
window  reflection  is  shown  in  Fig.  6(b).  The  relative  output 
phase  ripple  of  ±0.1  radians  is  consistent  with  (11).  However, 
in  the  case  when  the  current  source  is  dependent  on  the  cavity 
RF  fields,  the  electron  bunching  current  can  undergo  a  phase 
shift  and  interacts  with  the  reflected  wave.  In  this  case,  the 
maximum  magnitude  of  the  relative  output  phase  can  be 
estimated  from  (5)  to  be 

\Aip 

max  l^Atan  (rr^)-  (12> 

In  other  words,  an  output  phase  ripple  of  up  to  ±0.2  radians  can 
be  expected  for  a  -20-dB  reflection. 

Thus  far,  MAGY  simulation  results  have  only  been  shown 
for  the  constant  RF  current  excitation  of  the  output  cavity.  To 
illustrate  the  influence  of  window  parameters  on  the  perfor¬ 
mance  of  gyroklystrons,  we  have  also  performed  MAGY  sim¬ 
ulations  of  the  four-cavity  W-band  gyroklystron  with  electron 
beams.  The  main  difference,  with  respect  to  the  constant  current 
source  cases,  are  that  electron  beams  enter  the  output  cavity  pre¬ 
bunched  under  the  influence  of  the  previous  three  cavities  and 
that  the  bunched  electron  current  can  be  further  modulated  by 
the  output  cavity  RF  fields.  Shown  in  Fig.  7(a)  is  the  MAGY 
simulated  band  widths  plot  for  a  6- A,  65 -kV  electron  beam  at 
saturation  with  and  without  reflection.  More  details  about  the 
device  will  be  discussed  in  later  paragraphs  and  can  also  be 
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Fig.  6.  (a)  Output  phase  versus  frequency  with  and  without  window  reflection  and  (b)  corresponding  relative  phase  ripple,  Aip,  for  constant  RF  current  source 
excitation  of  output  cavity  [Fig.  3]. 


found  in  [3]  and  [4].  For  these  simulations,  the  beam  pitch  angle 
(a  =  v±/vz)  is  1.22  with  no  velocity  spread  and  the  magnetic 
field  is  36.4  kG  with  a  —  1%  taper  over  the  circuit.  Note  that  the 
center  frequency  shiwn  in  this  plot  is  now  93.84  GHz,  reflecting 
an  upshift  of  220  MHz  from  the  cold  resonance  frequency  of 
93.62  GHz  [Fig.  3],  because  of  the  effect  of  beam  loading.  As 
previously  mentioned,  the  window  is  ~73.6  cm  from  the  output 
cavity  with  a  -20-dB  reflection.  The  power  ripple  is  plotted  in 
Fig.  7(b).  In  contrast  to  Fig.  3(b),  when  the  output  cavity  is  ex¬ 
cited  with  a  constant  RF  current  source,  the  ripple  magnitude 
here  is  less  than  20%  and  is  much  larger  at  the  band  edges  than  at 
the  center.  This  effect  has  been  previously  explained  in  the  dis¬ 
cussion  following  (7).  The  output  phase  for  this  case  is  plotted 


in  Fig.  8(a),  and  the  relative  output  phase  ripple  is  shown  in 
Fig.  8(b).  As  expected,  the  maximum  magnitude  of  the  phase 
ripple  is  larger  here  than  that  shown  in  Fig.  6(b),  but  is  less 
than  0.2  radians,  as  can  be  expected  from  (12).  Note  that  the 
frequency  spacing  between  ripples  is  the  same  as  in  the  con¬ 
stant  current  excitation  case,  as  expected. 

The  dependence  of  the  bandwidth  ripple  on  the  length  be¬ 
tween  the  output  cavity  and  the  window  is  illustrated  by  Fig.  9 
for  two  different  lengths  L  -  40  and  130  cm.  The  simulation  pa¬ 
rameters  are  exactly  those  described  in  the  previous  paragraph. 
Comparing  with  Fig.  7(a),  it  can  be  seen  that  the  frequency 
spacing  is  larger  (smaller)  for  shorter  (longer)  lengths  in  accor¬ 
dance  with  (8).  In  addition,  the  frequency  at  which  the  peaks 


Fig.  7.  (a)  Saturated  bandwidth  with  and  without  window  reflection  and  (b)  corresponding  power  ripple,  A P/P0,  for  a  6  A,  65  kV,  a  -  1.22  electron  beam. 
Beam  enters  output  cavity  prebunched  due  to  interactions  in  the  three  previous  cavities. 


and  valleys  of  the  bandwidth  ripple  occur  is  also  extremely  sen¬ 
sitive  to  the  position  of  the  window.  This  is  demonstrated  in 
Fig.  10  for  two  window  positions,  one-quarter  of  an  axial  wave¬ 
length  apart  (~0.8  mm).  The  interchange  between  peaks  and 
valleys  can  clearly  be  seen.  This  behavior  is  apparent  upon  in¬ 
spection  of  (4)  and  (6).  Dependence  of  bandwidth  ripple  on  the 
reflection  levels  is  shown  in  Fig.  1 1  for  — 20-dB  [Fig.  6]  and 
— 26-dB  reflections.  As  expected,  the  ripple  magnitude  is  lower 
for  smaller  reflections. 

In  order  to  confirm  that  window  reflections  are  indeed  the 
cause  of  the  experimentally  observed  bandwidth  ripples,  as  il¬ 
lustrated  in  Fig.  1 ,  a  series  of  MAG  Y  simulations  have  been  per¬ 


formed  for  direct  comparison  with  the  data  using  known  exper¬ 
imental  parameters  from  the  four-cavity  W-band  gyroklystron 
experiments  [3],  [4].  The  resonance  frequency  and  quality  factor 
for  each  of  the  four  cavities  in  the  simulations  are  tuned  within 
the  experimental  uncertainty  of  the  cold-test  measurement  [11]. 
The  simulated  geometry  of  this  amplifier  is  illustrated  in  Fig.  2. 
The  beam  pitch  angles  and  velocity  spreads  used  in  the  sim¬ 
ulations  are  extrapolated  from  EGUN  [12]  simulations  of  the 
MIG  [10].  The  extrapolation  is  based  on  experimental  start-os¬ 
cillation  data  and  MAGY-predicted  pitch  angle  for  the  output 
cavity  oscillation  threshold,  and  on  the  relative  change  of  gun 
voltages  and  magnetic  field  settings  used  in  these  experiments. 
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Fig.  8.  (a)  Output  phase  versus  frequency  with  and  without  window  reflection  and  (b)  corresponding  relative  phase  ripple,  Atp,  for  a  6  A,  65  kV,  a  =  1 .22  electron 
beam  [Fig.  7]. 


The  gun  and  interaction  circuit  magnetic  field  are  calculated 
with  the  magnetics  code  EFFT  [13]  using  the  experimental  cur¬ 
rent  settings  and  th£  superconducting  coil  dimensions  as  input. 
Measured  and  calculated  magnetic  field  profiles  are  in  excellent 
agreement.  Most  importantly,  a  window  located  ~73.6  cm  (ac¬ 
cording  to  the  design)  from  the  output  cavity  with  a  reflection  of 
—20  dB  (window  cold-test  value)  is  assumed  in  all  simulations. 

Shown  in  Figs.  12  and  13  are  comparisons  of  MAGY  sim¬ 
ulated  and  experimentally  measured  saturated  bandwidths  and 
drive  curves,  respectively,  at  four  different  magnetic  fields  for 
4-A,  65-kV  electron  beams.  The  agreement  between  measured 
data  and  simulated  results  is  very  good.  The  beam  pitch  angle 


and  transverse  velocity  spread  used  for  each  case  is  as  shown 
in  the  plots,  as  are  the  input  power  for  the  bandwidth  and  the 
drive  frequency  for  the  drive  curves.  Note  that  the  input  power 
given  here  is  from  the  measured  power  at  the  input  cavity  flange, 
which  includes  the  ~3-dB  loss  between  this  flange  and  the  TEoi 
cavity.  The  axial  magnetic  field  profile  employed  in  the  simu¬ 
lations  is  shown  in  Fig.  14,  as  calculated  using  EFFT.  Similar 
simulation/measurement  agreement  is  also  obtained  for  band¬ 
width  and  drive  curves  with  6-A,  65-kV  electron  beams  at  four 
magnetic  field  profiles,  as  illustrated  by  Figs.  15  and  16,  respec¬ 
tively.  The  corresponding  magnetic  field  profiles  are  as  shown 
in  Fig.  17. 
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Fig.  9.  Saturated  bandwidth  plots  with  and  without  reflection  of  a  6  A,  65  kV>  a  =  1 ,22  electron  beam  as  in  Fig.  7  but  for  two  different  distances  between  window 
and  output  cavity  (a)  L  =  40  cm  and  (b)  L  =  130  cm. 


The  good  agreement  between  simulation  and  the  experiment 
demonstrated  here5|  together  with  the  absence  of  bandwidth 
ripple  in  simulations  with  no  reflection,  strongly  suggests  that 
window  reflections  are  indeed  the  cause  of  the  observed  band¬ 
width  ripple  in  the  high-average  power  W-band  gyroklystron. 
In  addition  to  the  bandwidth  ripple,  window  reflections  also 
induce  output  phase  ripple,  as  discussed  earlier.  Even  though 
the  discussion  given  here  has  been  restricted  to  window  reflec¬ 
tion,  the  analysis  is  general  and  can  be  applied  to  other  types  of 
reflectors.  Moreover,  we  have  shown  that,  similar  to  the  large 
signal  code  MAGYKL  [14],  MAGY  can  be  used  for  designing 
gyroklystrons.  In  particular,  MAGY  application  is  required  for 


describing  self-consistent  effects.  This  is  the  topic  of  the  next 
two  sections. 


III.  Higher  Order  Mode  Excitation  in  Output  Taper 

High-average  power  gyroamplifiers  often  employ  an  output 
taper  to  transition  from  the  near-cutoff  dimensions  of  the  inter¬ 
acting  circuit  to  a  larger  dimension  waveguide  [15].  This  transi¬ 
tion  is  necessitated  by  thermal  considerations  of  the  spent-beam 
collectors  or  by  requirements  dictated  by  the  output  windows. 
Nonlinear  uptapers  are  commonly  preferred  over  discrete  steps 
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Fig.  10.  Sensitivity  of  bandwidth  ripple  to  window  positions.  The  two  rippled  curves  have  window  position  shifted  by  a  quarter  of  an  axial  wavelength  (A L  = 
0.8  mm).  A  window  reflection  of  —20  dB  and  a  6  A,  65  kV,  a  =  1.22  electron  beam  were  simulated  as  in  Fig.  7. 


Frequency  (GHz) 


Fig.  11.  Dependence  of  saturated  bandwidth  on  level  of  window  reflection  for  a  6  A,  65  kV,  a  =  1.22  electron  beam. 

power  W-band  gyroklystron  discussed  in  the  previous  section 
[Fig.  2].  This  uptaper,  radially  shaped  to  follow  a  Dolph-Cheby- 
shev  profile,  serves  as  the  transition  from  the  radius  of  the  TEoi 
output  cavity  (~0.2  cm)  to  the  radius  of  the  beam  collector 
(~1.62  cm)  in  a  length  of  approximately  10  cm.  The  beam  col¬ 
lector  is  also  the  output  RF  waveguide  [10]. 

A  more  detailed  picture  of  the  operation  of  this  uptaper  can  be 
seen  in  Fig.  1 8,  which  shows  the  geometry  and  the  RF  mode  con¬ 
tents  as  a  function  of  axial  length  from  aMAGY  simulation.  The 
simulation  was  performed  by  exciting  the  output  cavity  with  a 


or  linear  tapers  because  of  minimal  mode  conversions,  thus,  en¬ 
suring  high  mode  purity.  High  mode  purity  is  critical  to  min¬ 
imize  undesirable  RF  absorption  and  reflection  in  the  ensuing 
transmission  line  and  to  reduce  unwanted  side  lobes  of  radiated 
pulses,  when  coupled  to  an  antenna  system. 

Because  of  physical  length  constraints  in  most  devices,  all  up- 
tapers  are  carefully  designed  to  maximize  the  mode  purity  with 
minimal  length.  The  final  optimized  geometric  profile  of  the  up¬ 
taper  is  selected  to  satisfy  this  requirement.  An  example  of  such 
an  optimized  design  is  the  nonlinear  uptaper  in  the  high-average 
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Fig.  12.  Comparison  between  experimental  bandwidth  (dots)  and  MAGY  results  (solid  lines)  for  4-A,  65-kV  beams  at  four  different  magnetic  field  settings. 
The  corresponding  circuit  magnetic  field  profiles  are  as  shown  in  Fig.  14.  The  input  power  and  beam  velocity  ratio  and  perpendicular  spread  for  each  case  are, 
respectively,  (a)  170  W,  1.5,  4.3%,  (b)  90  W,  1.55, 4.0%,  (c)  250  W,  1.42, 4.6%,  and  (d)  340  W,  1.32,  5%. 

constant  RF  current  source,  which  has  an  axial  Gaussian  profile.  A  potential  major  perturbation  to  this  delicate  balance  is 
Five  different  modes  are  included  in  the  simulation,  i.e.,  TEqi,  the  interaction  of  higher  order  modes  with  the  electron  beam 
TE02>  •  ■  ■>  and  TE05.  As  discussed  previously,  this  is  equivalent  in  the  nonlinear  uptaper.  The  interaction  is  expected  to  be 
to  performing  a  “numerical”  cold  test  of  the  output  cavity  and  that  of  the  gyro-TWT  type.  The  seed  that  drives  these  higher 
the  output  taper.  Thp  resulting  final  higher  order  mode  content  order  mode  interactions  will  be  the  power  converted  from  the 
of  less  than  —37  dB  from  the  simulation  confirms  the  validity  TEqi -mode  during  the  transition  process.  For  each  mode,  this 
of  the  design.  However,  it  is  interesting  \o  note  the  presence  of  interaction  takes  place  approximately  where  the  gyro-TWT 
substantial  higher  order  mode  within  the  transition  [16].  For  in-  resonance  condition  is  met,  i.e.  the  intersection  between 
stance,  the  TE02  content  can  be  as  high  as  20%  before  being  the  Doppler-shifted  beam  line,  ,  and  waveguide  dispersion 

reconverted  back  to  the  TEqi -mode.  This  process  of  converting  line,  u>  =  £lc(z)  +  kz(z)vz ,  and  waveguide  dispersion  line, 

the  operating  mode  to  higher  order  modes  and  then  reconverting  u2  —  u)^(z)2  -f-  kz(z)2c2.  Here,  £lc  is  the  electron  rela- 

back  to  the  operating  mode  to  obtain  high  mode  purity  in  a  rea-  tivistic  cyclotron  frequency,  is  the  cutoff  frequency  of  the 

sonable  length  is  a  delicate  one.  Consequently,  it  can  be  ex-  TEon-mode,  and  kz  is  the  axial  wavenumber.  Fig.  19  shows 
pected  that  perturbations  to  the  RF  flow  within  the  uptaper  can  the  cutoff  frequency  of  several  modes  as  a  function  of  axial 
quickly  compromise  the  mode  purity  of  the  RF  output.  distance  for  this  particular  geometry.  Also  shown  in  the  figure 


t 
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Fig.  13.  Comparison  between  experimental  drive  curves  (dots)  and  MAGY  results  (solid  lines)  for  4- A,  65-kV  beams  at  four  different  magnetic  field  settings. 
The  corresponding  circuit  magnetic  field  profiles  are  as  shown  in  Fig.  14.  Drive  frequency  is  93.8  GHz  for  all  cases.  The  beam  velocity  ratio  and  perpendicular 
spread  for  each  case  are,  respectively,  (a)  1.5,  4.3%,  (b)  1.55,  4:0%,  (c)  1.42,  4.6%,  and  (d)  1.32,  5%. 


Fig.  14.  Circuit  magnetic  field  profiles  employed  in  MAGY  simulations 
shown  in  Figs.  12  and  13.  Cavity  positions  and  dimensions  are  also  indicated 
for  clarity.  Field  profiles  are  calculated  with  magnetic  field  code  EFFI  using 
experimental  current  settings  for  the  super-conducting  coils  as  input. 


is  the  line  for  94  GHz  for  reference  purposes.  Assuming  that 
the  operating  frequency  is  94  GHz,  it  can  be  seen  that  higher 
order  mode  gyro-TWT  interactions  can  only  take  place  inside 
the  uptaper  approximately  when  the  cutoff  frequency  falls 
below  the  operating  frequency.  Because  of  the  rapidly  changing 
radius  of  the  taper,  it  is  also  expected  that  the  interaction  for 
each  mode  will  only  occur  over  a  fairly  short  distance  from 
where  the  resonance  condition  is  first  approximately  satisfied. 

This  phenomenon  is  most  pronounced  in  regimes  in  which  the 
TEqi  interaction  in  the  output  cavity  is  not  very  efficient,  such  as 
in  the  linear  regime  (low  drive  power)  or  at  the  bandwidth  edges. 
This  is  because  in  such  cases  thebeam  is  still  energetic  as  it  exits  the 
output  cavity  and  can  strongly  interact  with  higher  order  modes  in 
the  output  taper.  To  quantify  this  phenomenon,  Fig.  20  shows  the 
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Fig.  15.  Comparison  between  experimental  bandwidth  (dots)  and  MAGY  results  (solid  lines)  for  6- A,  65-kV  beams  at  four  different  magnetic  field  settings.  The 
corresponding  circuit  magnetic  field  profiles  are  as  shown  in  Fig.  17.  The  beam  velocity  ratio  and  perpendicular  spread  for  each  case  are,  respectively,  (a)  1.17, 
5%,  (b)  1.24,  4.0%,  (c)  1.21,  4.5%,  and  (d)  1.22,  4.3%.  Input  power  is  340  W  for  all  cases. 


total  output  power  as  well  as  the  power  distribution  in  higher  order 
modes  as  a  function  of  axial  distance  with  low  drive  power  (~0.5 
W  at  93.85  GHz)  from  aMAGY  simulation.  Thebeam  currentand 
voltage  for  this  simulation  are  6  A  and  65  kV,  respectively,  with  a 
pitch  angle  of  1 .3  and  no  velocity  spread.  A  uniform  magnetic  field 
of  36  kG  has  also  been  assumed.  The  resulting  TEoi  -mode  power 
from  the  output  cavity  is  about  1.8  kW,  and  it  can  be  readily  seen 
that  the  higher  order  modes  are  dominant.  The  total  output  power 
is  6.5  kW.  It  can  also  be  seen  that  the  lower  the  radial  mode  num¬ 
bers,  the  closer  the  interaction  to  the  output  cavity,  as  expected. 
Furthermore,  the  higher  order  mode  powers  also  peak  quickly, 
indicating  the  region  of  synchronism  for  each  mode  is  short.  A 
similar  simulation  with  zero  input  power  gives  zero  output  power. 
Thus,  confirming  that  higher  order  mode  excitation  is  indeed 


driven  by  the  seed  provided  by  the  mode  conversion  process  in  the 
output  taper. 

TE0i  and  higher  order  mode  powers  as  a  function  of 
frequency  are  plotted  in  Fig.  21(a)  at  the  saturation  of  the 
TEoi-mode  (drive  power  is  380  W).  It  can  be  seen  that  in  this 
case  the  higher  order  mode  content  is  much  reduced  relative  to 
the  total  output  power.  This  is  illustrated  in  Fig.  21(b),  where 
the  ratio  of  higher  order  mode  to  total  output  power  is  plotted 
as  a  function  of  frequency.  This  is  because  the  TEoi  interaction 
is  efficient  in  robbing  the  beam  of  much  of  its  transverse  kinetic 
energy.  Consequently,  higher  order  mode  interactions,  which 
occur  later,  are  much  less  efficient.  The  parameters  for  this 
particular  simulation  are  as  in  the  previous  case  but  with  a  = 
1.22  and  a  realistic  magnetic  field. 
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Fig.  16.  Comparison  between  experimental  drive  curves  (dots)  and  MAGY  results  (solid  lines)  for  6- A,  65-kV  beams  at  four  different  magnetic  field  settings. 
The  corresponding  circuit  magnetic  field  profiles  are  as  shown  in  Fig.  17.  The  beam  velocity  ratio  and  perpendicular  spread  for  each  case  are,  respectively,  (a) 
1.17,  5%,  (b)  1.24,  4.0%,  (c)  1.21,  4.5%,  and  (d)  1.22,  4.3%.  Drive  frequency  is  93.8  GHz  for  all  cases. 
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Fig.  17.  Circuit  magnetic  field  profiles  employed  in  MAGY  simulations 
shown  in  Figs.  15  and  16.  Cavity  positions  and  dimensions  are  also  indicated 
for  clarity.  Field  profiles  are  calculated  with  magnetic  field  code  EFFI  using 
experimental  current  settings  for  the  super-conducting  coils  as  input. 
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Fig.  1 8.  Geometry  of  output  cavity  and  nonlinear  uptaper  for  the  high  average 
power  W-band  gyroklystron.  Also  shown  are  axial  power  profiles  of  the  TE0i, 
TE02 ,  and  TE03  modes,  as  well  as  the  total  power  from  all  modes  with  no  beam. 


This  study  indicates  that  the  beam-wave  interaction  must  also  purity;  i.e.,  other  taper  profiles  can  reduce  the  amplification  of 

be  taken  into  account  in  uptaper  designs  to  ensure  high  mode  the  internal  higher  order  modes.  In  addition,  rapid  tapering  of 
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Fig.  19.  Geometry  and  cutoff  frequency  versus  axial  distance.  Also  shown  is  the  94  GHz  line  for  reference. 
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Fig.  20.  Mode  contents  of  the  output  power  versus  axial  distance  in  uptaper  at  low  drive  power  (0.5  W  at  93.85  GHz)  with  beam-wave  interactions  in  uptaper. 
Beam  current,  voltage,  and  pitch  angle  are  6-A,  65-kV,  and  1.3,  respectively. 


the  magnetic  field  in  the  uptaper  region  can  also  reduce  these 
interactions.  Experimental  confirmation  of  this  phenomenon  is 
planned. 

IV.  Beam-Induced  RF  Excitation  in  Cutoff  Drift 
Section 

Gyroklystrons  employ  drift  sections  in  part  to  RF  isolate  the 
cavities.  These  drift  sections  are  completely  cutoff  to  the  oper¬ 
ating  mode  at  the  operating  frequency  band,  and  they  are  usually 
lined  with  lossy  dielectric  to  suppress  spurious  oscillations.  For 


instance,  in  the  high-average  power  TEqi  W-band  gyroklystron, 
lossy  dielectric  linings  on  the  drift  tube  wall  are  employed  to 
suppress  potential  oscillation  of  modes  with  field  distributions 
similar  to  the  TEu-mode  in  circular  waveguide.  In  this  case,  the 
dielectric  causes  the  field  profiles  to  become  hybrid  (i.e.,  they 
have  both  axial  magnetic  and  axial  electric  field  components). 
However,  with  the  appropriate  choice  of  dielectric  properties 
and  radial  thickness,  the  hybrid  HEn-,  HE12-,  and  EHn-modes 
can  be  selectively  loaded  without  seriously  disturbing  the  cutoff 
properties  of  higher  order  modes  and  their  radial  field  profiles 


882 


IEEE  TRANSACTIONS  ON  PLASMA  SCIENCE,  VOL.  28,  NO.  3,  JUNE  2000 


Frequency  (GHz) 

(b) 

Fig.  21 .  (a)  TE0i  (left  scale)  and  higher  order  mode  (right  scale)  powers  as  a  function  of  frequency  at  saturation  and  (b)  ratio  of  higher  order  mode  power  to  total 
output  (in  dB)  for  a  6-A,  65-kV,  a  =  1.22  beam.  A  reflection  of  —  20-dB  has  been  assumed  for  all  modes. 


(within  the  vacuum  region  of  the  drift  tube).  For  a  complete  de¬ 
scription  of  this  procedure,  the  readers  are  referred  to  other  pub¬ 
lications  [4],  [17].  Because  the  objective  of  the  lossy  dielectric 
is  to  suppress  the  incipient  oscillation,  it  is  commonly  expected 
that  the  RF  fields  will  be  virtually  nonexistent,  and  thus,  thermal 
issues  will  not  be  a  serious  concern  even  during  average  power 
operation  in  the  design  mode.  Self-consistent  simulations  with 
MAGY,  however,  reveal  that  substantial  RF  levels  can  in  fact 
be  induced  in  these  drift  sections.  Interestingly,  the  induced  RF 


fields  are  found  to  be  that  of  the  operating  mode  and  are  at  the 
drive  frequency,  which  is  well  below  cutoff.  Moreover,  the  in¬ 
duced  RF  amplitude  is  found  to  scale  linearly  with  the  electron 
bunching  current. 

An  example  of  this  phenomenon  has  been  alluded  to  earlier 
in  connection  with  Fig.  2.  It  can  be  seen  there  that  below  cutoff 
excitation  is  evident  in  the  drift  section  between  the  penultimate 
cavity  and  the  output  cavity,  where  by  design  the  bunching  cur¬ 
rent  is  near  maximum  for  efficient  RF  extraction  in  the  output 
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cavity.  It  is  worth  mentioning  that  the  TE0i  cutoff  frequency  for  ( 
equivalent  all-metal  drift  sections  shown  is  approximately  118 
GHz,  and  that  the  operating  frequency  (hence,  the  induced  RF 
frequency)  is  about  93.8  GHz.  Note  that  in  the  absence  of  the 
bunched  beam,  the  RF  fields  for  the  operating  mode  with  just 
the  cavity  excited  are  evanescent  in  these  sections,  as  expected. 
However,  the  presence  of  the  bunched  electron  beam  can  excite 
RF  in  a  process  analogous  to  that  of  a  forced  harmonic  oscil¬ 
lator.  The  correspondence  can  be  explained  from  the  RF  wave 
equation 


In  MAGY,  at  the  present  time,  the  expected  level  of  thermal 
loading  to  the  dielectric  layers  caused  by  nonresonant  excitation 
is  estimated  from  the  radial  component  of  the  Poynting  vector 
evaluated  at  the  inner  radius  of  the  dielectric,  Rd-  For  TE-modes, 
this  is  given  by 

Pr(u,)=^Re{EeB:}r=Rd 

=  {Zs\Bz  l}r=-Rci  (16) 

where  Zs(r  =  Rd)  =  Eq/Bz  is  the  TE-modes  surface 
impedance  for  the  dielectric  layer.  In  the  case  when  the  outer 
radius  of  the  layer  is  also  the  conducting  wall  radius,  by 
matching  the  boundary  conditions,  the  TE  surface  impedance 
in  the  long  wavelength  approximation  can  be  shown  to  be 


1  d2E(t,r,oj) 


at 2 

47 r  dJb(t,r,w) 
c2  dt 


-V2£(i,r»  +  2 


47r<75  dE{t,r,ui) 


dt 


(13) 


where  E  is  the  induced  electric  field  in  the  drift  section,  Jb  is 
the  electron  bunching  current,  w  is  the  drive  angular  frequency, 
and  <js  is  the  surface  conductivity.  Using  the  long  axial  wave¬ 
length  approximation  (k2  =  0)  and  defining  em  as  the  normal¬ 
ized  transverse  eigenvector  for  the  mode  of  interest  (e.g.,  TEoi), 
the  amplitude  of  the  electric  field,  Em,  is  given  approximately 
by 


$ Eyn{t ,  Z j  Cj)  UJC  dErn(t)Z)U}') 


dt2 


+ 


Qs 


dt 


+  u)cErn(t,  z^uf) 


J_  Hl(l)n(xd)  +  pTEHf\xd) 
^  y/H&\xd)  +  pTV,Hn\xd)  . 


with 


Pte  =  -H'(1\xw)/H'(2\xw) 


—U) 

Xd.v:  —’\/£—Rd,u 


(17) 


=  ,SmdA.  (,4) 

Here,  ujc  is  the  cutoff  frequency  of  the  drift  section,  Q$ 
represents  surface  loss,  and  the  integral  is  over  the  cross  section 
of  the  drift  section.  Equation  (14)  has  the  form  of  a  damped 
forced  harmonic  oscillator  (see,  for  example,  [18]),  for  which 
the  steady-state  solution  is 


and  e  is  the  complex  dielectric  constant  of  the  layer.  In  these 
equations,  and  are  Hankel  functions,  n  is  the  az¬ 
imuthal  mode  index  number,  and  the  prime  denotes  the  deriva¬ 
tive.  By  relating  the  induced  axial  RF  magnetic  field  Bz  to  the 
bunching  current  via  the  wave  equation 

Vj _Bz{r,u)  +  ~Bz{r,oj)  =  ±(Jb(r,u)xez)  (18) 


the  power  loading  onto  the  dielectric  layer  from  the  induced  RF 
-4niu  I  Jb(z:  Lj)  •  em  dA  fields  can  be  calculated  according  to  (16). 

Eni(z,cj)  — _  _  (15)  3"o  confirm  the  validity  of  this  phenomenon,  we  have  per- 

cj2  —  u)2  iiouc/Qs  formed  simulations  with  the  PIC  code  MAGIC  [19].  These 

simulations  also  reveal  the  existence  of  the  below  cutoff  exci- 
It  can  be  readily  seen  here  that  the  amplitude  of  the  induced  tation  phenomenon  in  the  presence  of  bunched  electron  beam, 
fields  scales  linearly  with  the  electron  bunching  current,  which  in  agreement  with  MAGY.  Additional  MAGIC  simulations 
explains  the  observable  amplitude  of  the  fields  in  the  drift  sec-  were  performed  to  estimate  the  level  of  power  loading  on 
tion  between  the  penultimate  cavity  and  the  output  cavity.  the  dielectric  layer.  These  MAGIC  results  have  been  com- 

Two  implications  are  associated  with  the  nonresonant  RF  ex-  pared  with  MAGY  estimates  using  the  formulation  of  (16) 
citation  at  the  drive  frequency  phenomenon.  The  first  is  that,  and  found  to  be  in  good  agreement.  Consequently,  (16)  is 
strictly  speaking,  the  usual  assumption  of  negligible  RF  fields  presently  used  for  estimating  the  thermal  impact  on  ceramics 
in  cutoff  sections  or  severs  in  amplifiers,  and  thus,  that  electrons  because  of  this  effect.  A  more  self-consistent  formulation  of 
only  undergo  ballistic  bunching,  is  not  correct.  As  shown  here,  lossy  materials,  however,  is  being  implemented  in  MAGY 
the  RF  field  amplitude  can  indeed  be  observable  depending  on  [20]. 

the  bunching  process  and  how  far  the  drive  frequency  is  from  To  quantify  the  effect  of  dielectric  power  loading  from 
cutoff.  More  importantly,  for  high-average  power  devices,  the  nonresonant  excitation,  we  have  applied  MAGY  to  one  of  the 
induced  RF  fields  in  the  drift  sections  can  thermally  load  lossy  W-band  gyroklystrons  at  NRL  [21].  The  geometry  and  RF 
dielectric  layers  employed  there  for  stability  purposes.  This  im-  field  profile  (at  saturation)  are  as  shown  in  Fig.  22.  In  fact, 
plies  that  careful  attention  must  be  paid  to  the  cooling  of  these  this  gyroklystron  is  the  prototype  for  the  high-average  power 
dielectric  layers.  gyroklystron  [3],  [4].  This  prototype,  as  in  the  high-average 
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high-average  power  device  has  been  performed.  The  analysis 
Fig.  23.  Comparison  of  operating  mode  radial  electric  field  profiles  in  drift  indicates  that  a  moderate  temperature  rise  of  25  °C  relative  to 
sections  with  and  without  the  dielectric  linine  at  93.8  GHz.  The  metal  wall  radius  ,  ,  •  ,  .  .  , 

is  0.175  cm.  The  dielectric  (a  composite  of  80%  BeO  and  20%  SiC  by  weight)  the  cooled  metal  wal1  can  be  achieved  assuming  good  thermal 

inner  radius  and  thickness  are  0. 1 75  and  0.08  cm,  respectively.  Cutoff  frequency  contact  between  the  ceramics  and  the  cooled  metal  wall, 

is  0.3%  lower  with  the  dielectric  lining. 


power  version  [4],  employs  dielectric  loading  in  each  of  the 
three  drift  sections  to  suppress  spurious  oscillations.  The 
dielectrics  are  a  composite  of  80%  BeO  and  20%  SiC  by  weight 
with  a  complex  dielectric  constant  e  =  12.24  -  i3.67  (at 

94  GHz).  The  dielectric  inner  radius  and  thickness  are  0.175 
and  0.08  cm,  respectively.  The  lengths  are  0.3  cm,  except  for 
the  one  in  between  the  penultimate  and  output  cavities,  which 
is  0.48  cm  long.  With  this  configuration  of  dielectric  loading 
in  the  drift  tube,  there  is  minimal  effect  on  the  field  profile 
and  cutoff  behavior  of  the  operating  mode,  which  is  similar 
to  the  TE0i  in  the  vacuum  region  and  has  a  small  tail  that 
leaks  into  the  dielectric,  as  illustrated  in  Fig.  23.  At  93.8  GHz, 
the  radial  shift  in  the  vacuum  field  profile  maximum  is  only 
about  0.3%,  with  an  equal  percentage  change  in  the  effective 
wall  radius  (resulting  in  a  cutoff  frequency  of  104  GHz).  This 
is  negligible,  justifying  the  use  of  the  metal  drift  tube  with 
unchanged  dimensions  and  an  effective  finite  conductivity.  The 
peak  power  loading  on  the  ceramics  as  calculated  by  MAGY 
for  a  6  A,  65  kV,  a  =  1.5  at  saturation  is  shown  in  Fig.  24. 
As  expected,  the  loading  is  highest  for  the  last  dielectric  layer 
because  the  electron  bunching  current  is  near  maximum  there. 
The  induced  RF  field  at  this  position  is  indeed  observable,  as 
can  be  seen  in  Fig.  22. 

For  high-average  power  operation,  the  highest  peak  power 
loading  of  2.4  kW/cm2  [Fig.  24]  translates  into  an  average 
power  loading  as  high  as  264  W/cm2,  assuming  a  duty  cycle 


V.  Conclusion 

In  this  paper,  the  self-consistent  capability  of  MAGY 
is  being  employed  to  investigate  several  issues  relating  to 
gyroklystrons.  The  effects  of  window  reflections  on  the 
performance  of  gyroklystrons  were  presented.  It  was  shown 
that  window  reflection  could  induce  bandwidth  power  ripple. 
The  amplitude  of  which  depends  strongly  on  the  reflection 
coefficient,  and  the  ripple  frequency  spacing  is  a  function  of  the 
geometry  of  the  waveguide  between  the  output  cavity  and  the 
reflector.  Furthermore,  window  reflections  also  result  in  output 
phase  ripple  across  the  band  and  modify  the  output  cavity 
diffractive  quality  factor  Quiff  -  This  modification  of  Qu^  can 
significantly  impact  the  start  oscillation  current  threshold  of  the 
output  cavity.  The  excellent  correlation  between  experimental 
data  and  simulation  results  indicates  that  the  observed  ripple  in 
the  bandwidth  of  the  high-average  power  W-band  gyroklystron 
[3],  [4]  can  indeed  be  attributed  to  a  small  -20-dB  reflection 
from  the  window.  Simple  scaling  laws  are  also  developed, 
which  allow  for  a  quick  estimation  of  the  impact  on  the 
gyroklystron  performance  for  a  given  window  reflection.  The 
potential  for  higher  order  mode  excitations  in  nonlinear  output 
tapers  is  also  discussed.  The  higher  order  mode  interactions 
are  that  of  the  gyro-TWT  type  with  the  initial  drive  power  for 
each  mode  provided  by  the  internal  mode  conversion  in  the 
taper.  This  effect  can  compromise  mode  purity  of  the  output 
RF  for  which  the  use  of  output  tapers  is  intended  to  minimize. 
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It  is  most  pronounced  in  at  low  drive  power  or  at  the  bandwidth 
edges  where  the  operating  mode  interaction  in  the  output  cavity 
is  not  very  efficient.  This  suggests  that  designs  of  output  taper 
should  take  into  account  possible  beam-wave  interactions 
to  ensure  that  high  mode  purity  is  indeed  achieved.  Finally, 
nonresonance  excitation  in  cutoff  drift  sections  driven  by  the 
bunched  electron  beam  was  also  presented.  The  induced  RF 
fields  are  found  to  be  that  of  the  operating  mode  and  are  at  the 
drive  frequency,  which  is  well  below  cutoff.  The  RF  amplitude 
scales  linearly  with  that  of  the  electron  bunching  current.  This 
phenomenon  is  similar  to  that  of  a  forced  harmonic  oscillator. 
The  presence  of  induced  RF  in  drift  sections  can  potentially 
impact  the  use  and  design  of  lossy  ceramics  employed  in  these 
sections  for  stability  purposes,  particularly,  in  high-average 
power  devices  where  thermal  issues  are  a  major  concern. 

Even  though  the  work  presented  here  focus  exclusively  on 
gyroklystrons,  MAGY  has  also  been  applied  to  gyro-TWT  de¬ 
signs.  The  results  of  gyro-TWT  modeling  will  be  presented  in  a 
future  publication. 
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